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CHAPTER I 
INTRODUCTION 

The search for natural resources, the need for accurate weather 
prediction, as well as many other scientific studies have stimulated 
widespread interest in remote sensing devices. One such device is the 
microwave radiometer. The radiometer has proven valuable in remotely 
determining the physical temperature of surfaces, thickness of sea and 
lake ice, salinity of seawater, and the location of oil spills. The 
radiometer may be able to detemine these physical properties by mea- 
suring the power level of the electromagnetic energy emitted from the 
object. The energy emitted by the object is actually a very low level 
gaussian random variable with zero mean [i]-]. Estimating the variance of 
this signal presents several interesting problems. One problem is that 
the actual power level of the received signal may be lower than the 
receiver noise power. Because of the extremely low level of these 
signals, substantial amplification is required before any attempt can be 
made to measure the power level of the received signal. The amplifiers 
used to provide this gain must be very stable since changes in the gain 
will cause error in the measured power of the received signal. Thus, 
in order to measure the power of the received signal accurately, the 
receiver noise and the gain of the radiometer must be precisely known 
and extremely stable. Any uncertainty in the receiver noise power or 
the gain results in erroneous measurement of the power of the signal and 
thus will degrade the quality of the measurement of physical properties. 
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The Dicke switching noise injection radiometer shown in figure 1-1 
is one system used to estimate the power level of the received signal 
while reducing errors due to the previously mentioned problems. An 
introduction to the basic operation of this radiometer can be found in 
Appendix A and also in references [ 4 ], [ 5 ], and [7]. Referring to fig- 
■ure 1-1 j the received signal is modulated by switching between a refer- 
ence noise source and the antenna output. This Dicke switching strategy 
eliminates the effect of the receiver noise on the measurement of the 
received signal. The noise injection or feedback scheme eliminates the 
effect of gain fluctuation on this measurement. A Dicke switching noise 
injection radiometer of this type has been developed by R. F. Harrington 
at NASA Langley Research Center. Since construction of the radiometer, 
advances in understanding and interpretation of radiometric data have 
necessitated further refinements in the existing system. 

The objective of this research is to investigate techniques for improv- 
ing the power measurement performance of the existing Dicke switching noise 
injection radiometer. The general approach of this investigation will be to 
develop an analytical model of the existing radiometer. It will be 
shown from this model that the use of sample data techniques leads to an 
improved power estimation strategy. 

A model of the existing radiometer will be developed and analyzed 
in Chapter II. From this analysis features which adversely affect the 
power estimate performance of the radiometer will be identified. In the 
first section of Chapter III, this model is used to develop a sample data 
processing approach which will improve the power estimate performance 
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of the radiometer. The second section of Chapter III deals with the 
details of implementing the design improvements. The results of an 
analysis using the developed model and the proposed processing tech- 
nique, as well as preliminary experimental results obtained via a prototype 
design using an Intel 8086 microprocessor, are presented in Chapter IV. 
Conclusions drawn from these results are contained in Chapter V. 



5 


CHAPTER II 

DEVELOPMENT OF RADIOMETER MODEL 
In this chapter an- analytical model will be developed for the 
existing Dicke switching noise injection radiometer. Several performance 
parameters of the radiometer, such as sensitivity and measurement time, 
will be defined. From the analytical model, expressions for these per- 
formance parameters will be obtained. Methods of improving the power 
measurement performance of the radiometer will then be discussed. 

A block diagram of the existing radiometer is shown in figure 2-1. 
The system consists of five major sections: the antenna subsystem, the 

microwave front end, the detector, the video amplifier, and the analog 
signal processor. Of interest here is the development of a model to be 
used in subsequent analysis of the entire radiometer. In order to 
evolve such a model, the transfer functions of the various sections must 
be determined. To provide continuity of presentation, several brief 
descriptions of the subsystems are presented. A detailed analysis of 
these subsystems has been performed and can be found in reference [lO]. 

In the following development of the radiometer model, the signal and 
noise power components will be represented by their respective effective 
temperatures. The basis for this representation is developed in 
Appendix A. Henceforth, the term effective temperature will be used to 
refer to a corresponding power. The brightness temperature of an object 
will be used to refer to the temperature corresponding to the power 
received by the antenna. 
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Figure 2-1, Block diagram, of existing radiometer 
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A. Development of AnaXytical Model 

Referring to figure 2-1, the antenna subsystem collects the inci- 
dent electromagnetic radiation and feeds it to the microwave front end. 
The feedback noise is also injected into the microwave front end and is 
added to the received signal. The circulator alternately switches the 
input of the tunnel diode amplifier between the received signal and the 
injected noise combination and the reference noise. The output of the 
tunnel diode amplifier is down converted to zero intermediate frequency 
(i.e., dc homodyne receiver) and amplified further by the mixer pre- 
amplifier. This entire microwave section of the radiometer is housed 
in a constant temperature enclosure. The reference noise is provided by 
a matched load also within the enclosure. Thus, the physical tempera- 
ture of the microwave components and the reference temperature are both 
at the enclosure temperature of 308 K. This eliminates the need to 
accurately measure the ohmic losses within the microwave front end. The 
frequency response of the microwave front end is sufficiently broad that 
the effect of this portion of the system can be represented by a constant 
power density gain factor. This gain bandwidth constant was determined 
to be 137.63 dB-Hz [lO]. 

The output of the microwave front end is filtered and detected by 
the detector section. The filter bank of figure 2-1 establishes the 
effective input statistical bandwidth which will be of importance in 
later analysis. Since changing this input filter will clearly change 
the power applied to the square-law detector, attenuators are included 
in the filter bank to assure that the power applied to the detector is 
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the same for all filters. The output of the square-law diode is a 
linear function of the input power. Measurements indicated that the 
gain constant of this square-law detector is itOO V/W [lO]. 

The relatively low level square— law detector output is then ampli- 
fied hy the video amplifier. This amplifier shoiild pass virtually all 
components of the Dicke switch frequency with minimum attenuation and 
phase shift. The amplifier need not pass dc since the mean of the 
temperature estimate appears as a modulated component of the Dicke fre- 
quency. Thus, the amplifier should display substantial attenuation at 
low frequencies to reduce flicker or 1 /f noise. This video amplifier 
has been analyzed [lO], and the transfer function was found to be 

G (3) = 0.693(1 + 2 . 39 s)s^ 

^ [1 + s/(2tt X 2.lt-l^[^l + s/(2tt X 2.0l)J 

X —i ( 2 “l) 

[1 + s/( 2 tt X o. 3 )J [1 + s/( 27T X 5.03 X io 3 )j 

An analysis of the Bode plot of this function revealed that the video 
amplifier is relatively flat from above a few hertz to about 5 or 
6 kHz. 

The output of the video amplifier is applied to the correlator, 
which effectively demodulates the signal restoring the dc value. The 
output of the correlator is then applied to the analog signal processor 
or estimation and averaging section. The analog processor must estimate 
the mean or dc level of the output of the demodulator switch. This 



9 


estimation section consists of an integrator and a loop filter. The 
transfer function for this combination was determined to be 


1'96.4[i + s/(2tt X 62.529^ 

s[l + s/(2tt X 5.684^ 


The final constant required to form the complete transfer function 
of the radiometer is the feedback constant determined by the noise 
injection process. The noise injection is accomplished by the use of a 
voltage-tO“frequency converter producing constant width pulses with a 
frequency directly proportional to the output voltage of the analog 
processor. The output pulses of the voltage-to- frequency converter con- 
trol a PIN [13] diode switch which is used as a gate controlling the 
noise diode. The noise diode excess noise temperature, or noise above 
ambient temperature, is specified by the manufacturer to be 31 dB above 
290 K, corresponding to an excess noise temperature of approximately 
365,000 K. The voltage-to-frequency converter produces an output of 
863 pulses/second per volt. Each pulse is TO ys wide; thus, the duty 
cycle is given by 


d = (863) (TO X 10-^) Vq (2-3) 

d = O.060U Vq 

where is the voltage output of the estimation and averaging section. 

These injected noise pulses are attenuated by a 6 dB attenuator and then 
added to the antenna noise via a 20 dB directional coupler. The noise 
diode and attenuator are both included in the constant temperature 
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enclosure with the microwave front end. Thus, they have a physical 
temperature equal to the reference temperature, i.e., 308 K. Consider 
the effect of the 6 dB attenuator on the feedback noise temperature. 

The output of an attenuator at a physical temperature Tp with input 
noise temperature Tjjj can be shown to be 

% " “ L)Tp {2-k) 

where L is the power loss factor (L < l). The input to the attenuator 
is the total noise diode temperature Tp = Tpj^ + Tp. Therefore, since 
the attenuator is within the constant temperature enclosure, equa- 
tion (2-i+) can be written as 

Tp = 0.25Tp + (1 - 0.25)Tp 

'^N ~ ^ ■*■(!” 0.25)Tp (2-5) 

\ = 0.25TJ3, + T^. 

Thus, an eo,uival6n‘t modGl for "the feodhack can be obtained by including 
only the excess noise temperature of the noise diode and adding the 
reference temperature at the output of the attenuator. Using this 
method, the feedback mechanism shovn in figure 2-2 is obtained. 

The feedback model of figure 2—2 can be combined with the system 
constants calculated earlier to form the block diagram model of the 
radiometer shown in figure 2-3. The input bandwidth determined by the 
filter bank of figure 2-1 need not be considered in the control loop 
analysis. However, its effect on the operation of the radiometer will 
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Figure 2-2. Illustration of noise addition model. 
























be considered later. Note also that the frequency response of the video 
amplifier is assumed to he sufficiently broad not to affect the control 
loop, and it is represented by its pass band gain. 

The control loop of figure 2-3 can be simplified by first combining 
the feedback constants, then moving Boltzmann's constant within the con- 
trol loop, and finally combining the gain constants in the forward loop. 

A further simplification is obtained by representing the switching pro- 
cess with only one switch. This is possible since fluctuations in the 
frequency response of the system between the first switch and the corre- 
lator switch occur at such high frequencies that they can be neglected 
for control loop analysis. With these simplifications, the control model 
illustrated in figure 2~k is obtained, where the 0.5 factor is due to 
the Dicke switching process. 

Before analyzing this model, an important point about the feedback 
pulses and the video amplifier should be noted. The assumption was made 
that the frequency response of the video amplifier was sufficiently 
broad to be neglected in the control loop analysis. However, the ampli- 
fier is not broadband enough to pass the narrow feedback pulses without 
some distortion. Proper operation of the Dicke switching process 
requires that all noise injected into the antenna portion of the Dicke 
cycle appear in the antenna portion at the correlator switch. Any energy 
that "spills over" due to bandwidth limitation between the switches will 
result in an error in the measurement process. It was determined that 
the bandwidth of the video amplifier was insufficient to process the 
noise injection pulses without some "spillover" error. This error was 
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observed qualitatively using computer simulation of the radiometer. The 
error has also been measured during calibration of the radiometer [lO]. 

An analysis of the control model of figure 2-k will now be conducted 
in order to develop a closed-loop transfer function for the present 
radiometer. The noise temperature Tg of figure 2-4 is given by 

T^ = 0.5(Ti - Tfi) (2-6) 

where Tj is the effective output temperature of the directional coupler 
and T-g is the reference temperature. The effective output temperature 
of the directional coupler can be seen from figure 2-4 to be 


Tj = (1 - + Kjcdg T„) (2-7) 

where is the injected noise appearing at the attenuator output. 

The noise temperature at the attenuator output is given by 


% %Tt'^EX%^o 


(2-8) 


where is the attenuator constant, Tg^^ is the excess temperature 

of the noise diode, and is the constant relating the output 

voltage to the duty cycle of the injected noise pulses. The output 
voltage Vq(s) can be expressed as 

Vo(s) = -G(s)Tp. (2-9) 

Substituting equations (2-6), (2-T), and (2-8) into (2-9) results in the 
following expression for Y^{s)t 
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V^(s) = - I G(s)(Tj - Tg) 

I G(s)[(l - Kdc)Ta + %c(Tb + %) - Tg] 

( 2 - 10 ) 

I G(s) [(1 - Kbc)Ta + Kbc(Tb + KattTexKiVq) - TbJ 
- G(s)|jl - Kbc)(T^ - Tb) + KBcKATT'^EX%''^q] 

where the term ^dc^ATt'^EX^I^o feedback noise. The transfer 

function relating Vq(s) to Tb - T^ can be determined from equa- 
tion (2-10) to be 


V^(s) = - 
Vo(s) = - 
V^(s) = - 


Gcl(s) - 


Vq(s) 

^B-^A 


i G(s)(1 - Ktsp) 


r (-) - ci ■ DC 

Gqb(s) 

1 + ^ G(s)KB(,K^rp^Tg^Kj^ 


( 2 - 11 ) 


Substituting the proper values into (2-ll) results in 


|(0.99)(1.132) 

1 + 2.5l^5 X I0"^s 


s(l + 28 X l0“3s) 

1 + |(0.99)(1.132) 

1 + 2.51+5 X I0"^s 

55.12 

_s(l + 28 X io"3s)_ 


(2-12a) 


^ , s 20.0(1 + 2.5^5 X 10-3s) 

(j-prls; - - ^ 

3^=^ + 38.52s + 1103.1 


(2-1 2b) 


0 


0 




The poles of this transfer function are readily determined to be 
s = -19.26 ± J27.06. Thus, the damping factor is a = 19.26, and the 



damped radian frequency is = 27.06 rad/s. The undamped natural 
radian frequency is Wq = 33.21 rad/s. The damping ratio 
C = 0.580. This closed-loop transfer function can he used to 

analyze several aspects of the radiometer's performance. 

B. Analysis of Model 

A thorough understanding of exactly what the radiometer is measuring 
and how this measurement relates to the brightness temperature are essen- 
tial in attempting to properly • interpret the results obtained using the 
previously developed model. As explained in Appendix A, the output of 
the square-law detector is a random process with mean proportional to 
the brightness temperature of the scene within the antenna footprint. 

Let y(t) be this process, and let 

y(t) = y (t) + z(t) (2t13) 

where z(t) is a random process with zero mean, and y(t) is the mean 
value of y(t) corresponding to the brightness temperature within the 
antenna footprint. The brightness temperature of the scene B(t) is 
then determined by 

B(t) = E{y(t)} (2-lU) 

where B(t) is the brightness temperature of the scene within the 
antenna footprint versus time. Since it is not practical to have a 
large number of radiometers measuring the same scene in order to form 
the ensemble average of equation (2-li|), the radiometer must use a time 
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average of the output of the square-law detector as an estimator of the 
brightness temperature. It is shown in Appendix B that this time 
average is not an unbiased estimator of the brightness temperature, but 
it is an estimate of the average brightness temperature of the scene 
covered by the path of the antenna footprint in T seconds. Thus, 
longer integration times will produce a more accurate estimate of the 
average brightness temperature of a larger area. If, however, more 
spatial resolution of the actual brightness temperature is required, 
longer integration times may not provide more accuracy. If the time 
average is used as an estimator of B(t), and if the integration time 
is comparable to the antenna footprint dwell time, a bias error results. 
This error was calculated in equation (B— it) and is repeated here for 
convenience. 


OO 



n=l 


(2n) 
B(t) T 


2n 


(2n) !(2n + 


( 2 - 15 ) 


This bias error or "antenna smear" error is due to the antenna footprint 
movement and subsequent changing of the brightness temperature during the 
integration time. It should be noted that the bias error of equa- 
tion (2-15) is a function of the brightness temperature. Thus, the 
level of the "antenna smear" error will depend on the type of scene 
being measured. 

As the analog processor forms the required time average and pro- 
duces an estimate of the brightness temperature, the random nature of 
the process causes fluctuations in the estimate. The sensitivity of a 




a 



19 


radiometer is an indicator of the level of these fluctuations. The 
sensitivity AT of a radiometer is defined as the amount by which the 
input temperature must be changed to result in the output estimate 
changing by one standard deviation. The sensitivity of a closed-loop 
square wave correlated Dicke radiometer can be shown to be [12], 


AT = 2(Tj^ + Tg) 



(2-l6) 


where Tg and Tg are the receiver and reference noise temperatures, 
respectively, is the effective input statistical bandwidth, and 

%0 output equivalent noise bandwidth. Thus, in order to deter- 

mine the sensitivity of the radiometer, the output equivalent noise 
bandwidth must be known. 

The equivalent one sided noise bandwidth of the radiometer is 
defined by 

|GgL(f)|^df (2-17) 



where GQg(f) is the response of the control model given in equa- 
tion (2-12). It can be shown [lO] that the equivalent one sided noise 
bandwidth for the transfer function of (2-12) is given by 


1 + 


8C 




%0 = 


(2-18) 
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where T2 is the time constant of the lead portion of the loop filter 
and is the undamped natural radian frequency. Substituting the 

appropriate values from equation (2-12) into equation (2-l8) results in 


B. 


1 + (33.21 X 2 . 5 h 5 X lp-3)^ 


NO 

Bjjq = 7 • 21 Hz . 


8(0.580) 


33.21 


(2-19) 


Substituting this value for Bj^q as well as 308 and 627 K for Tg 
and Tg, respectively, into equation (2-l6), the following result is 
obtained: 


( 2 - 20 ) 

AT = 1870 

V ®SI 

The various sensitivities as a function of different input handwidths 
were calculated using equation (2-20), and the values for the four 
available system input bandwidths are summarized in table 2-1. 


AT = 2(627 + 308) , - 

V ®SI 


TABLE 2-1 

SENSITIVITY AND REQUIRED MEASUREMENT TIME 
FOR PRESENT RADIOMETER 


%I 


20 MHz 

100 MHz 

500 MHz 

2 GHz 

AT 


1.59 K 

0.710 K 

0.317 K 

0.159 K 

^S 


0.393 s 

0.409 S 

0.419 s 

0.511 S 




There is a settling time associated with the loop parameters used 
to obtain the output noise bandwidth of equation (2-19). In order to 
determine this settling time or the time required to form an estimate, 
the transient response of the loop to a step input must be obtained. 

The step response of the present radiometer can be obtained by multi- 
plying equation (2-12b) by the Laplace transform of a, step function and 
taking the inverse transform. The step response y(t) is determined 
to be 

y(t) = C(T£ - T^)[l - l.lTe-19-26t (27.o6t + 0.956)] 

Clearly, the error Ay(t) due to the loop settling time is 

Ay(t) = C(Tg - T^)(l.l7)e--'^9.26t g^^^ (27.06t + O.956) 

If the settling time is defined as the time required for the settling 
error to be less than 10 percent of the fluctuation of the estimate AT, 
then equation (2-22) can be written as 

II = (Tb - TA)(l.l7)e-19.26t g.^ (27.06t + 0.956) (2-23) 

where AT is the sensitivity of the measurement and is the time 

required to make the measurement. When determining the maxim\jm settling 
time from equation (2—23), the peaks due to the overshoot phenomena must 
be considered. The time at yhich these peaks occur can be found by 


( 2 - 21 ) 


( 2 - 22 ) 



taking the derivative of equation (2-23) and setting it equal to zero. 
This yields 

(2-2k) 

The settling time was calculated for each input bandwidth using equa- 
tions (2-23) and (2-24), and these results are presented in table 2-1, 
Increasing this settling time will decrease the fluctuations AT, but 
will also increase the possibility of "antenna smear" error. Thus, for 
optimum performance the radiometer should form the estimate to the 
required sensitivity in the shortest possible integration time. 

In addition to the noise fluctuations and "antenna smear" error in 
the output of the loop, there is an error component due to the Dicke 
switching process. This so-called Dicke ripple is due to the modulation 
and demodulation processes within the Dicke radiometer. Consider the 
block diagram of figure 2-5. If it is assumed that the loop has reached 
steady state and- if the feedback pulses are represented by their average 
level, the voltage waveform of 2-5 is obtained for the output of the 
correlator. Since this waveform, has odd symmetry, it has no even har- 
monics and the magnitudes of the odd ntunbered Fourier components are 
readily determined to be 

% = + Tp)|G(jN0J^)| 


max 


27.06 


TT 


(2n - 1) 


- 0.956 


) 


n — 1,2,3, 


(2-25) 







2k 


where to^ is the radian Dicke frequency. The dc level Vq with the 
loop in steady state is given by 


= 


Ko 


( 2 - 26 ) 


For typical loop parameters, the fundamental component V2_ is usually 
much larger than the higher order components. Thus, the assumption that 
the total Dicke ripple level can be approximated by will be made. 

This assumption was verified through computer simulation. The ratio of 
the fundamental component of the Dicke ripple to the dc level at the 
output of the processor can then be determined to be 



P = 


4 ^r) 

TT (Tb - Tp) 


K^KgG(jW^) 


(2-27) 


Substituting the appropriate value for the loop gain into equation (2-27) 
results in 


8 (308 + 627 ) 

^ TT (308 - T^) 

= 3.658 X 100 ^ 

308 - 

This output peak ripple level expressed as a percentage of the output 
signal is shown as a function of antenna signal temperature in fig- 
ure 2-6. From this curve, it is readily observed that the ripple level 


31.195(1 + 2.545 X 10 “-^Jw^) 


Jw^d + 28 X 10 ~ 2 jw^) 


X 100^ 






f 





increases markedly with increasing input signal temperature and becomes 
quite significant for large antenna signal temperatures. 

This Dicke ripple appears about the dc level at the output of the 
analog processor. Thus, as the ripple level increases the dynamic range 
available for the dc level decreases. This is undesirable since the dc 
level is the actual estimate of the received power. As the ripple level 
Increases the dynamic range required of the Injection noise source 
Increases. If the ripple becomes excessive the temperature range over 
which the radiometer can operate may be reduced. 

Through the use of the system model and the transfer function of 
equation (2-12), several aspects of the radiometer’s performance were 
investigated above. The sensitivity of a radiometer was defined, and 
an expression for the sensitivity of this particular radiometer was 
obtained. The ripple associated with the Dicke switching process was 
calculated. The implications of time averaging a nonergodic process 
were examined. An expression for the bias error resulting from the use 
of a time average as an indicator of the ensemble average of a nonergodic 
process was developed. It was noted that this bias or "antenna smear" 
error became larger for longer averaging times, and that the bias error 
was a function of the actual brightness temperature to be measured. 

Thus, if the time average is to be used as an estimator of the actual 
brightness temperature, the optimm averaging or integration time will 
be a tradeoff between sensitivity and "antenna smear" error. It was 
determined that for a specified spatial resolution the optimum integra- 
tion time is a function of the brightness temperature of the scene 



bsin^ iH6asured, Thus, th© ©sisimatiion sch©m© can be improved by perform- 
ing measurements of several integration time and sensitivities. How- 
ever, care must be taken since as the loop is made to respond faster 
to reduce the "antenna smear" error of a quickly changing scene, the 
Dicke ripple may increase, driving the radiometer into saturation. 

Thus, the power estimation strategy can be improved by reducing the 
Dicke ripple within the control loop and developing a variable integra- 
tion time averaging scheme. In the following chapters a digital pro- 
cessing unit is designed, developed, and implemented which will improve 
the performance of the radiometer. 
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CHAPTER III 

PROPOSED PROCESSING STRATEGY ^ 

In this chapter a power estimation strategy is developed which will 

1 * 

improve the tradeoff between sensitivity and "antenna smear" error. It 
will- be shown in this chapter that digital processing techniques can be 
used to improve the power estimate performance of the radiometer. 

Anomalies of the sample data approach will be exploited to reduce the 
Dicke ripple which will aid in the development of a variable integration 
time post-loop processing strategy. This post-loop strategy coupled 
with the added flexibility of the digital signal processor will result 
in an improved power estimation strategy. 

A. Improved Processing Technique 

When applying digital signal processing techniques to the radiom- 
eter, one of the first questions to be answered is at what point within 
the radiometer should the analog signal be sampled. Sampling as close 
to the antenna as possible seems desirable since a larger number of 
analog components will be replaced by the digital processor which is 
more stable and flexible. However, sampling prior to the square-law 
detector is not possible with presently available components since system 
bandwidths of 1 GHz or greater are required for proper radiometer opera- * 

tion. System bandwidth between the square— law detector and correlator 
should be maintained above 5 kHz for proper operation. However, after 
the square-law detector only the mean of the signal is of interest, and 
the question arises as to whether the signal can be undersampled without 



undesirable effects. It has been shown, however, that although recon- 
struction of the signal in the conventional sense is not required, 
sampling below the Nyquist rate will cause a noise foldover effect [lO]. 
This foldover effect is the conventional aliasing error, and it reduces 
the mean to variance power ratio. Therefore, if the signal is sampled 
before the correlator, a sample rate above 10 kHz would be required. It 
may be possible to develop the required digital processor with a data 
rate of 10 kHz, thus placing the aliasing filter and analog-to-digital 
(A/D) converter before the correlator. However, after the correlator, 
the signal could be filtered any amount since at this poipt any low-pass 
filtering provides additional smoothing in the estimate of the dc level 
of the signal. This fact coupled with some technical difficulties of a 
digital processor operating at 10 kHz resiilt in the conclusion that 
sampling after the correlator is most desirable. 

The selection of the sampling frequency should be made very care- 
fully. Due to the periodic nat\ire of the Dicke switch, the output of 
the correlator has line components at harmonics of the Dicke switch fre- 
quency. If the sampling frequency is selected as an odd harmonic of the 
Dicke switch frequency, the line components will be translated to dc by 
the sampling process causing a substantial error in the measureijient of 
the actual dc level. However, due to the symmetry of the deterministic 
signal at the output of the correlator (with no injected noise), no line 
components exist at even harmonics of the switch frequency. Depending 
on the feedback method, the feedback noise pulses "riding" on this square 
wave may cause line components to be formed at these even harmonics. 
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Even so, sampling at an even harmonic of the Dicke switch frequency may 
still he very desirable. 

Let us assume that an accumulator will exist within the micropro- 
cessor algorithm to replace the analog integrator. Consider the dis- 
crete transfer function R(z) of the accumulator. 




= -J cot (^) ( 3 _ 2 ) 

■where is the sample interval. The amplitude response R(e'^ 

is shown in figure 3-1. This amplitude response approximates that of an 
ideal integrator for f « 1/2 Tq. Note, however, that for f = 1/2 t 

s 

the steady-state response is identically zero. That is, the accumulator 
has a transmission null at one-half the sampling frequency. Herein lies 
the motivation for selecting the sampling frequency as a harmonic of the 
Dicke switch frequency. If the sampling frequency is chosen as twice 
the Dicke switch frequency, the transmission nulls of the accumulator 




will completely eliminate the line components at all odd harmonics of 
the Dicke switch frequency. Thus, all components of the Dicke ripple, 
except those introduced during the feedback process, have been elimi- 
nated. This should greatly reduce the undesirable Dicke ripple. It is 
also worth noting that if the sampling frequency is twice the Dicke 
frequency, the Dicke switch and analog-to-digital converter can easily 
be synchronized and can be controlled by the same clock. 

If the sampling frequency is to be twice the Dicke frequency, an 
optimum Dicke frequency must be determined. The selection of the Dicke 
frequency is a compromise between ease of aliasing filtering and an 
excessive bandwidth requirement for amplifiers in the forward loop of 
the radiometer. A higher Dicke frequency provides more guard band 
between dc and harmonics of the switching frequency, thus easing aliasing 
filter requirements. However, a higher Dicke switch frequency also 
requires wider bandwidth amplifiers in the forward loop since these 
amplifiers should pass all significant components of the Dicke frequency. 
Through computer simulation, a Dicke switch frequency of 300 Hz proved 
to provide acceptable aliasing filter requirements and did not require 
excessive bandwidths of the forward loop amplifiers. If the Dicke 
switch frequency is chosen as 300 Hz, the sampling frequency should be 
chosen as 600 Hz. The accumulator then eliminates odd harmonics of the 
Dicke switch frequency, thus reducing the Dicke ripple effect. 

The aliasing filter is selected to minimize the aliasing or noise 
foldover problem discussed earlier without substantially affecting the 
control loop; i.e., the closed-loop response should be exclusively 



determined by the digital processor algorithm. With the selected 
sampling frequency of 600 Hz, the aliasing filter should attenuate the 
spectrum above 300 Hz to minimize noise foldover. The only portion of 
the signal that is of interest at this point is dc, so the flatness of 
the passband is not of interest. A simple low— pass filter having three 
real poles at 100 Hz was found to meet these requirements. At 300 Hz, 
this network attenuates the spectrum 30 dB and rolls off at a rate of 
60 dB/decade. Computer simulations indicated that this network reduced 
the noise foldover to insignificant levels, and the effect of this 
aliasing filter on the closed-loop response was found to be minimal. 

It was noted earlier that line components in the correlator output 
may be formed at even harmonics of the Dicke switch frequency due to the 
feedback pulses. If such line components were formed, they would be 
translated directly to dc when sampled at 600 Hz and would cause sig- 
nificant error. Simulations proved that all these even harmonic com- 
ponents except the one at 600 Hz were attenuated sufficiently to cause 
negligible error. The component at 600 Hz is attenuated kj dB by the 
aliasipg filter. However, when sampled, it is translated to dc and 
causes an error in the measurement. To eliminate this error, a notch 
filter with a null at 600 Hz is added to the aliasing filter. 

The sampled output of the aliasing filter is operated on by the 
microprocessor algorithm. The algor j.thm must form an estimate of the 
mean of the input signal. Just as for the analog processor, the problem 
of finding an integration time for minim-um error arises. As discussed 
in Chapter II, a slower control loop, i.e. , increased integration time. 



will reduce fluctuation error, but will increase the “antenna smear" 
error or reduce spatial resolution* The digital processor, however, 
can adjust the loop time constant under program control. Thus, an 
operator, or possibly the microprocessor, itself, could estimate the 
"antenna smear" error given in equation (2-15) and appropriately adjust 
the integration time to minimize total error for a specified spatial 
resolution. To provide this flexibility, the algorithm must estimate 
the mean of the input over several integration times. There are two 
obvious methods which would provide this capability. First, the loop 
filter could be adjusted under program control, thus modifying the loop 
response. The second method would be to design a first order loop with 
a settling time equal to the shortest integration time required. Post- 
loop averaging can then be performed on the loop output to reduce esti- 
mate variance. It was decided that the second method would be most 
desirable since the first order loop has certain design advantages, and 
the estimates for all accuracies are available simultaneously. 

In order to form the required first order loop, the loop filter con- 
tained in the analog process is eliminated from the microprocessor algo- 
rithm, and the analog integrator is replaced by an accumulator. As 
mentioned earlier, the accmulator has the desirable effect of elimi- 
nating the odd harmonics of the Dicke ripple. The output of the accumu- 
lator will control the feedback noise, thus completing the loop. It was 
determined that an integration time of 100 ms was the minimum value 
required. A first order loop with a damping factor of 50 will provide 
this minimimi integration time. The loop gain required to obtain a 



35 




damping factor of 50 will be calculated in the next chapter. To provide 
the longer integration times, the output of the loop is further averaged 
by a post-loop processor. 

The post-loop processor sums N loop samples and then divides 
by N. To simplify the algorithm, N is restricted to integer powers 
of two. To limit the dynamic range req.uired and avoid the possible need 
for floating point arithmetic, a partitioning averaging scheme will be 
implemented. In this averaging scheme, which is actually several sum- 
and-dump algorithnjs , two samples out of the loop are added together and 
then divided by two. This output will be denoted as yg(n), where n 
is the post-loop sampling integer. Then y 2 (n) is added to the previous 
value, i.e., y 2 (n— 2), and divided by two resulting in yj^(n). Let 

y^(n) represent discrete values appearing at the output of the loop. 

The remaining values out of the post-loop processor can be expressed as 


y2(n) 

yij(n) 

y^Cn) 


y^(n) + y^(n - l) 

2 

yi(n) + y3^(n - l) + y^(n - 2) + y^(n - 3) 

- j- 

y2(n) + y2(n - 2) 

2 


= I ^ - i) 

i=0 

yi^(n) + yi^(n - 4) 

yg(n) = 


(3-3) 


(3-4) 


(3-5) 


2 
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yi6(n) 

yi6(n) 



73(11) + 73(11 - 8) 
2 


31 

i=0 


732(11) = 


yi6(n) + yi6(n - 1 ^) 
2 


63 

y6!+(^) = ^ 

i=0 

722 ^^^ '*' y 32 ^^ ~ 32) 
F- 


(3-6) 


(3-7) 


(3-8) 


The partitioning averaging scheme uses the previously averaged output to 
form the next averaged output. This partitioning scheme, in addition to 
reducing dynamic range requirements, uses far less memory and processor 
time than the direct averaging method. 

The post-loop averaging should decrease the AT of the loop output. 
To determine how well the sum-and-dump algorithm decreases the fluctua- 
tions of the measurement, the effect of this algorithm on the variance 
of the loop output should be studied. It will be assumed that the input 
to the post-loop processor is ergodic and stationary (negligible 
^antenna smear** error). The variance of the yj^(n) output of the post- 
loop processor is then determined to be 
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= Var I yjj( n) [ = Var 




(3-9) 


It can fee shown that if the samples y^^Cn - i) are statistically inde- 
pendent, equation (3-9) can fee written as 




( 3 - 10 ) 


Thus, if the samples at the loop output are statistically independent, 
the sum-and-dump averaging reduces the variance of the output fluctua- 
tions fey 1/N. 

In addition to providing data for the post-loop processor, the 
output of the first order loop must control the duty cycle of the 
injected noise pulses, thus controlling the feedback noise. Due to the 
sensitivity of the steady-state error of the radiometer to variations 
in the feedback path, it is desired to keep the feedback path within 
the processor algorithm as much as possible. The feedback will be 
accomplished fey loading the N bit word out of the loop into a down 
counter and clocking the down counter at a rate of 2-1 times the 
sample rate. The injected noise is on when the down counter contains a 
word greater than zero and is cutoff when the down counter reaches zero. 

In this way, the level of the injected noise power is adjusted to any 

N . 

of 2 levels. It is also desired to improve the injection scheme fey 
injecting the noise in the early portion of the Dicke cycle to reduce 
the loss of the injected p\ilse due to rounding of the pulses as dis- 
cussed in Chapter II. If the down counter is started at the beginning 
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of the antenna portion of the Dicke cycle, and if the maximum duty cycle 
is kept somewhat less than 100 percent, the injected noise loss due to 
rounding and spillover can be virtually eliminated. 

The maximum duty cycle or maximum injected noise corresponds to 
the minimum input temperature. Referring to the feedback path of fig- 
ure 2-2, based on a reference source of 308 K, a directional coupler 
coefficient of -20 dB, an excess noise diode temperature of 365,000 K, 
and the original 6 dB fixed attenuator, the maximum duty cycle is 

[(365,000) (0.25)] = 0.334 (3-11) 

In order to increase the length of a typical noise pulse and reduce 
quantization error in the feedback, the maximimi duty cycle is increased 
to near unity. To increase the duty cycle, the fixed attenuator was 
changed to 10 dB, and the maximum duty cycle then becomes 

= [o^“ 308j y^[( 365, 000) (0.10] = 0.835 (3-12) 

This noise injection scheme differs considerahly from the original design 
in that it will he injecting noise continuously during a portion of the 
Dicke cycle. The injection scheme provides excellent feedback stability ^ 
and with a maximum duty cycle of 0.835 will virtually eliminate the 
noise spillover into the reference half cycle. However, it does intro- 
duce a slight error due to the finite bandwidth of the video amplifier. 
The video amplifier has a bandwidth of 5 kHz and will reject the portion 
of the input spectrimi above 5 kHz. Since narrower pulses will have a 




wider frequency spectrum than wider pulses, the video amplifier will 
reject a larger proportion of power for narrow pulses than for wide 
pulses. This slight error can he compensated for by careful calibration 
and subsequent adjustment of recorded data. 

In summary, the proposed digital signal processor should filter the 
output of the correlator with three simple one pole sections, each with 
a cutoff frequency of 100 Hz, and a notch at 600 Hz. The ovtput of this 
aliasing filter is sampled at 600 Hz, and the samples are processed by 
the digital filter algorithm. The output of the digital filter algo- 
rithm is used to control the feedback pulses and is further processed by 
the post-loop processing algorithm. A block diagram of this proposed 
system is shown in figure 3-2. 

B . Implementation 

The hardware required to implement this proposed digital processor 
can be considered as being composed of four sections: an aliasing 

filter, an analog-to-digital converter, the microprocessor for imple- 
mentation of the algorithm, and a section to generate the timing signals 
required by the radiometer and to control the feedback. 

The aliasing filter is composed of the three low-pass sections and 
the notch filter. The three real poles of the aliasing filter ape pro- 
vided by the circuit shown in figure 3^3. The transfer function of this 
circuit can be shown to be 
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Figure 3-3- Three pole filter. 
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For the component values shown in figure 3-3, equation (3-13) hecomes 


GFj_(s) = 


(200ir)- 


[s^ + s(800ir) + (200tt)^ [s + (2007t)3 


(3-l4) 


Gp, ( s) - 


/l + — - — \ 
\ 200TT ) 


From equation (3-lU), the circuit of figure 3-3 clearly provides the 
required three real poles. The notch at 6 OO Hz is developed using the 
circuit illustrated in figiure 3-U. The transfer function associated with 
figure 3-4 can be shown to be 


GF2 (s) 


+ (I200ir)^ 

s^ + s(4800ir) + (I200ir)^ 


( 3 - 15 ) 


Combining the three pole and notch filters of figures 3-3 and 3-4 results 
in the required aliasing filter shown in figure 3-5, where IC2 is 
inserted to isolate R^C^ from the notch filter. When equations (3-l4) 
and ( 3 - 15 ) are combined the transfer function G(s) of the aliasing 
filter (shown in figure 3-5) is 


G(s) 


[s^ + (l200TT)^](200Tr)^ 

[s^ + s(4800it) + (12007t)2 (s + 200tt)3 






>: 


(3-16) 
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Figure 3-5. Required aliasing filter. 
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The aliasing filter requirement specified earlier in this chapter is 
met by equation (3-l6). Thus, the network shown in figure 3-5 is an 
adequate aliasing filter. 

The output of this filter must now be sampled by the analog-to- 
digital converter. The first aspect of this analog-to-digital converter 
to be considered is that resolution uncertainties should be significantly 
less than the fluctuations inherent in the measurement. Using equa- 
tion (2-l6) and the equivalent output noise bandwidth of the loop given 
in equation ( 4 - 8 ), the fluctuation of the output of the loop is given by 

( 3 - 17 ) 

where is the input statistical bandwidth. Substituting the 

appropriate values for Bgj the fluctuation AT at the output of the 
loop ranges from 2.1 K to 0.21 K corresponding to a range of B„-. from 
20 MHz to 2 GHz. The minimum input temperature of interest is 100 K; 
therefore, the fluctuations in the measurement at the loop output range 
from 2.1^ to 0.21^ of the values to be measured. The resolution of an 
N bit analog to digital in a bipolar offset format is given by 

r = ±2"®^ X 100 % (3-l8) 

Due to the fluctuations inherent in the measurement process as well as 
uncertainties in the radiometer exterior to the signal processing sec- 
tion, a resolution of approximately 10^ of the minimum AT (i.e., 

0.02^) was determined to be tolerable. Using equation (3-l8), the 
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number of bits required for this resolution was determined to be 12,1, 
This value is so close to 12, it was determined that 12 bits would be 
sufficient. The Analog Device's AD572 was selected as the analog-to- 
digital converter since, in addition to its 12 bit resolution, it is 
extremely linear and has a fast conversion time (25 ys). 

In order to obtain an accuracy within one-half of the least sig- 
nificant bit (±1/2 LSB) from the AD5T2, the input signal must not change 
by more than 1/2 LSB during the conversion time. The maximum amount by 
which the input can change during the conversion interval is given by 



(3-19) 


where is the input range of the analog-to-digital converter. If 

the range of the AD572 is adjusted to ±2.5 V (Vr = 5 V) then since the 
AD572 is a 12 bit converter AV is determined to be 


AV 



= 0.6l mV 


( 3 - 20 ) 


Since the conversion time of the AD572 is approximately 25 ys, the 
maximum allowable rate of change of the input voltage is 

^ _ 0.6l mV 

max 25 ys 


V 






= 2h.k V/s 


(3-21) 



A more meaningful result may be obtained by computing the maximum fre- 
quency of a ±2.5 volt sine wave which does not exceed the maximum rate 
of change. The maximum rate of change of a sine wave occurs as it goes 
through zero and is equal to 


§1 

dt 


max 


(Vp)(2Trf) 


( 3 - 22 ) 


where Vp is the peak level of the sine wave. Solving for frequency 
in equation (3-22), and substituting 2.5 for Vp, the result is 


max 


2 _ 

5TT dt 


max 


= j^(2k.k) 


= 1.55 Hz. 


(3-23) 


Since the aliasing filter has a cutoff frequency well above 1.55 Hz, a 
sample-and-hold device must be included in the design. 

The requirements on the microprocessor are not very stringent since 
care was taken in developing the digital algorithm to minimize the 
arithmetic operations which the microprocessor would have to perform. 

The microprocessor must implement the accumulator and post-loop algo- 
rithms, and it will operate in fixed point arithmetic. The micropro- 
cessor will interpret the binary point as leading all bits, i.e., the 
value Is There Is a 1, to 1 expansion (two bits) 

through the accumulator, and the microprocessor must operate on l4 bit 
words to accommodate this expansion. Many microprocessors are available 
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which could easily fulfill these requirements with a data rate of 
600 Hz. An Intel 8086 microprocessor was used in the prototype design 
solely due to availability. The 8086 is a I6 bit microprocessor and so 
can easily operate on the required l4 bit words. However, with the 
relatively slow sample rate and simple algorithm, an 8 bit micropro- 
cessor operating with double precision could also implement the 
algorithm. 

To reduce the possibility of an overflow due to the expansion in 
the accumulator, the samples out of the analog-to-digital converter are 
sign extended and loaded into the 12 least significant bits of the 
microprocessor. Since the 8086 has 16 bit words,. the maximum input 
voltage of 2.5 V will be interpreted by the microprocessor as 2 ~^ - 2~^^ 
Thus, there is a "redefinition" gain of 


■'ado = (2-‘‘ - 2-^h/2.5 
= 0-025 


(3-24) 


The output of the acc\jmulator will be a l4 bit word located in the 
l4 least significant bits. In order to have available the most signifi- 
cant 12 bits, the output word is shifted right 2 bits which is equivalent 
to the divide by 4 shown in figure 3-2. When these 12 bits are sent to 
the down counter, it will interpret all 12 bits being "one" as the 
level l-2~^^. Thus there is another "redefinition" gain of I6. The 
control loop model including "redefinition" gains is shown in figure 3-6. 
A flow chart of the required algorithm is shown in figure 3-T- The 






Figure 3-6. Control model -of 


proposed digital processing radiometer 
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algorithm consists of two adds and two right shifts for the accumulator 
and six adds and right shifts for the post-loop processor, in a minimum 
of lit hit arithmetic. The microprocessor must also have available 
enough Read Only. Memory (ROM) to store the required instruction code and 
six l4 hit words of Random Access Memory (RAM) to store the previous 
values of the post-loop outputs. 

The digital processor must also generate the Dicke switch drive and 
the control pulse for the noise diode. These timing signals are obtained 
through the use of the Intel 8253 Programmable Interval Timer. The 8253 
consists of three independently programmable l6 bit counters. Each 
counter can be programmed to operate in any of six modes. Counter 0 of 
the 8253 will be programmed as a divide by N counter. In this mode, the 
output of counter 0 will stay high for N clock cycles, then go low for 
one clock cycle, and return high for another N cycles. This counter 
will be used to divide the 2.h^ MHz microprocessor clock bjy 2^^ - 1. The 
output of this counter will be a 598.29 Hz pulse train used to control 
the, analog-to-digital converter . Counter 1 will be programmed as a 
variable monostable multivibratpr (commonly referred to as a one-shot). 

In the one-shot mode, the output, of counter 1 will go low when the 
gate 1 input goes high and stays low for N clock cycles, where N is the 
digital word loaded into the counter. The output of the Ipop algorithm 
will be loaded into counter 1. It is then triggered at the beginning 
of the Dicke cycle. Thus, the duty cycle of the output of counter 1 is 
controlled by the output word from the accumulator algorithm. The third 



counter is also programmed as a one-shot and it is used to complete 
the "handshaking" between the anal6g-to-digital converter and the 
microprocessor. 

Combining the filter, analog-to-digital counter, and the 8253, the 
digital signal processor shown in figure 3-8 is formed. The output of 
counter 0 is inverted and used as the start convert pulse. The output 
of counter 0 is also used to clock the D-type flip-flop of the 
which converts the pulse train into a square wave of one-half the pulse 
train frequency. The square wave output of the flip-flop is used as the 
Dicke switch drive. As the analog-to-digital converter receives the 
start convert pulse and begins the conversion, the status line goes 
high. This status line enables the sample and hold circuit which holds 
the signal at the input to the analog-to-digital converter constant for 
the entire conversion interval. At the end of the conversion, the status 
line goes low which triggers the one-shot of counter 2 of the 8253. The 
output of counter 2 goes low directing the microprocessor to service the 
input port. The output of the loop algorithm is loaded into counter 1 
of the 8253 , which is programmed as a variable one-shot. Counter 1 is 
triggered by the Dicke drive such that the output goes low at the 
beginning of the antenna portion of the Dicke cycle and stays low until 
the counter content reaches zero. The output of counter 1 is used to 
control the noise diode, i.e., the injected noise. Note that all 
required timing is derived from the microprocessor clock (see figure 3-8) 
Thus, variations of the Clock frequency do not disturb the operation of 


the radiometer. 




Figure 3 8, Digital signal processor circuit diagram 
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The digital processor illustrated in figure 3-8 has heen developed 
with the idea of exploiting advantages of the digital process in mind. 
This processor should substantially reduce the Dicke ripple within the 
control loop, thus making possible integration times much shorter than 
were obtainable with the original design. The shorter integration times 
allow most of the averaging to be done in the more efficient post-loop 
algorithm. The flexibility of the digital processor and the use of the 
post-loop processing strategy produce simultaneous estimates of several 
integration times. This flexible averaging scheme allows an operator, 
or possibly the microprocessor itself, to select the integration time 
resulting in the best estimate for a specified spatial resolution. 
Results of an analysis of this new estimation strategy and experimental 
results obtained using the prototype discussed in this chapter are pre- 
sented in the next chapter to verify the improved estimation performance 


of the radiometer. 
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CHAPTER IV 
RESULTS 


An analysis of the proposed estimation strategy and the correspond- 
ing analytical results are presented in the first section of this 
chapter. Jn the second section experimental. /;j’esults will be presented. 

A. Analytical 


In order to determine a modei for the first order control loop, 
the discrete time transfer function of the accumulator given by 


R(z) 


1 + z 


- 1 / 


1 - z 


-1 




where z represents the Z-transform variable, must be transformed to 
its continuous time approximation. The direct relationship between the 
s-plane and the z-plane is given by 

z = e"'"" (M) 

where Tg is the sample interval. 

Then by use of the Taylor expansion, eq^uation (4—1) can be approxi- 
mated at low frequencies as 


R'(s) 


2 




(4-3) 


For the selected sampling rate of approximately 600 Hz, equation (4-3) 
becomes 


R'(s) 


1200 


s 


(4-4) 
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Referring to figure 3-6 and using the same values for the gain of the 
RF subsystems, since R'(s) is the dominant pole in the loop response, 
the forward loop transfer function is ^ 

(2.88 X 10-3)(l200)(KF)(0.1t)(lA) ^ 

G(s) = ^ 

where K^, is an adjustable gain in the forward loop. The feedback 
constant K 2 can be determined by considering how much noise is fed 
back for an output digital word of all "ones." When all "ones" appear 
at the output 5 the equivalent decimal value is interpreted to he 1-2"* 

The binary point is interpreted as leading all bits. When the output 
word contains all ”oneSj" the noise diode remains on for the entire half 
of the Dicke cycle, and the feedback noise is equivalent to 365 K. Thus, 
the feedback constant is K2 = 365. The closed-loop response; of the pre- 
posed digital processor as part of the radiometer can now be determined 
as 


Kp(l200)(2.88 X 10-3)(0.4)(l/h) 


Kp(l200)(2.88 X 10“3)(365)(0.U)(l/h) 

S 


Hcl(s) 


1 

1 + s . . 

126 . 1 k Kp 


(i+-6) 
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It was determined that a damping factor of 50 is required. Clearly, if 

a damping factor of 50 is desired, then should be 0.396, and (4-6) 

becomes 


HpL(s) 



( 4 - 7 ) 


The output of this loop is transferred to the post— loop processor. 
The post-loop processor then further averages the samples from the loop, 
and if the samples are statistically independent, it reduces the vari- 
ance of the mean of the samples by l/N. However, as the samples in the 
post-loop processor become statistically dependent, they provide less 
additional information about the process, and the variance of the esti- 
mate is reduced by less than l/W. 

To determine if the samples out of the loop are statistically 
dependent, the noise bandwidth of the loop should be found. The equiva- 
lent noise bandwidth of the output of the fir^t order loop is calculated 


as 



= 12.5 Hz 


( 4 - 8 ) 



Since the output equivalent noise bandwidth of the closed-loop response 
is 12.5 Hz, samples closer together than approximately l/(2Bjj]^) or 
ho ms apart will display a degree of correlation. The loop operates at 
600 Hz so samples are available at the output of the loop every 1.66 ms. 
Clearly, these samples would be highly correlated. Since the efficiency 
of the post-loop process decreases as the samples become correlated, it 
would be more efficient for the post-loop processor to sample the output 
of the loop such that it only operates on uncorrelated samples. Of 
interest here is the optimum reduction in variance in a specific averag- 
ing time. Thus, as the sample rate is reduced to insure uncorrelated 
samples, care must be taken not to reduce the number of available uncor- 
related samples such that the variance of the sample set is increased. 

/ When attempting to sample the output of the loop, the sample rate 
should be chosen to provide the maximum number of uncorrelated samples 
available from the loop output in the specified time. Oversampling, that 
is, taking samples too close together, will provide the best possible 
reduction in variance for the averaging time. However, this is not 
optimum since the same reduction in variance could have been obtained 
using fewer samples. Undersampling, that is, taking samples so far apart 
that fewer than the available uncorrelated samples are taken, will 
increase the variance of the sample set. Thus, although the samples 
are uncorrelated the increased variance of the sample set results in 
less than optimum reduction in the variance for the specified averaging 
time. It was decided that if the loop output were sampled at 37*5 Hz 
or every l6th sample out of the loop all available uncorrelated samples 
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would be obtained. However, at this sample rate (26.66 ms) the output 
will be oversampled and so the samples into the post-loop processor will 
be somewhat correlated. 

A model must now be developed describing the effect of the post- 
loop processor on samples whicl> are not statistically independent. In 
the earlier description of the post— loop processor, the assumption was 
made that the samples were statistically Independent. With this 
assumption the post-loop processor was found to reduce the variance of 
the samples by l/N. If this assumption is not made, the purely sta- 
tistical approach of equation (3-10) to modeling the post-loop processor 
becomes unwieldy, and the results obtained provide little insight into 
the process. 

An analysis of the post-loop processor from a frequency domain view- 
point will provide an adequate model of the post— loop processor. The 
post-loop processor is illustrated from the z-transform point of view in 
figure ^—1. Consider the form of the post— loop algorithm 

N-1 

V'') = I y yi(" - 0-9) 

^0 

where yjj(n) is the N averaged post-loop output and y^^(n) is tjie 
loop output. Let Yjj(z) represent the z-transform of y|j(n), and let 
Y^(z) represent the z-transform of y^^(n). Transformation of both 
sides of equation (U- 9) yields 

W-1 

' I 2 

i=0 


(U-10) 








6l 


The transfer function of the post-loop processor can be expressed as 


P(z) = 


Y N-1 

l(z) N Z 

i=l 


-i 




The steady-state frequency response can be found by substituting 

jtOTo 


z = e 


into equation (U-ll) 


N-1 


p(j») = I 2 




(4-12) 


i=l 


where is the post-loop sample time, which is equal to sixteen times 
the loop sample interval; i.e., = l6 T^. With the aid of the summa- 
tion formula 

N-1 . . N 


ST' i 1 - a 


a 


i=0 


(4-13) 


equation (U-12) can be written as 

JwTqN 


P(jO)) 


1 - e 


N(l - 


jWTo 


) 


P(f) = 


P(f) = 


1 - e 


J 2iTf t^N 


H(1 - 


-jTTfToN jiTfToN +jirfToN 
(e -_e )e 

^(e-jTrfXo _ ^jTTfTo^^+jTTfTo 


(4-14) 


+j(N-l)TTfTo 

TVT / -mS ^ 


N sin (iTT^P) 
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The eq^uivalen't one-sided noise bandwidth of the post— loop processor is 
given by 


B. 


NP 




1/2To 


|p(f) 1^ df 


(irNT^f ) ~]^ 
sin (TTTQf) 


^1/2Tq si 

"Jo F 


df 


(^-15) 


It can be shown [lO] that the value of the integral in equation (4-15) 
is given by 


B = _i_ 
NP 2 NTq 


(4-16) 


If a flat input noise power spectrum n is assumed, the input variance 
is given by 


a. = 
in IN 


n 

2^0 


(4-17) 


Substituting equation (4— l6) into equation (4-17) 3- result equivalent to 
the result of equation (3-10) is obtained. 



riB. 


NP 


n 

2Nt 




N 


(4-l8) 
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V 


w 




Clearly, the assumption of a flat input power spectrum at this point in 
the frequency domain approach is equivalent to the assumption of sta- 
tistical independence of samples in equation (3-10). 

The input post-loop samples are not statistically independent , or 
from the frequency domain point of view, the input power spectrum is not 
flat. The input noise spectrum is shaped before reaching the post-loop 
processor by the power transfer function of the loop. The loop relative 
amplitude-squared frequency response function A^(f ) has been determined 


to be 



(U-19) 


In order to determine the effect of the post-loop processor on sta- 
tistically dependent samples, the effect of the loop process A^(f) on 

Li 

the power spectrimi of the post-loop processor input must be examined. 

The noise bandwidth of the combination of the loop response and the post- 
loop processor is found to be 




fl/2To 

sin (NiTfTQ) 

'0 

N sin (TTftQ) 



(4-20) 


df 
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Letting x = T^f equation (4-20) becomes 



1 [ ~ sin (MTrx)~ j 

NTq Jq l_sin Ux) J 



(4-21) 


Equation (4-21) was evaluated numerically, and the results are shown in 
table 4-1 for selected values of N. 


TABLE 4-1 

EFFECTIVE NOISE BANDWIDTHS OF POST-LOOP PROCESSOR OUTPUTS 


N 


1 

2 

4 

8 

16 

32 

64 

^N 0 ^"o 


0.3326 

0.3459 

0.4177 

0.4586 

0.4793 

0.4896 

0.4948 

®N0 


12.5 Hz 

6.49 Hz 

3.92 Hz 

2.15 Hz 

1.12 Hz 

0.574 Hz 

0.290 Hz 


The sensitivity of the proposed digital signal processing radiom- 
eter can be determined using the expression for the sensitivity given 
in equation (2-l6), and the output equivalent noise bandwidths for the 
control-loop and post-loop processor are listed in table 4-1. Frqm 
equation (2-l6), the sepsitivity AT of the closed-loop noise injection 
radiometer is given as 

1 2B 

AT = 1870.-:^ (4-22) 

V ®SI 

where :/.s the output equivalent noise bandwidth and is the 

input statistical bandwidth. The sensitivity of the proposed design can 
now be found by substituting the output noise bandwidths for the various 




post-loop outputs, listed in table 4 - 1 , into equation ( 4 - 22 ). These 
sensitivities are listed in table 4 - 2 . 


^ TABLE 4-2 

SENSITIVITY AND REQUIRED MEASUREMENT TIME FOR 
DIGITAL PROCESSING RADIOMETER 


®SI 

N 

1 

2 

4 

8 

16 

32 

64 

20 MHz 

at (in K) 

2.09 

1.51 

1.17 

0.867 

0.626 

0.448 

0.318 

Tg (IN s) 

o.i6o 

0.187 

0.240 

0.320 

0.533 

0.960 

1.79 

100 MH^ 

AT (IN K) 

0.935 

0.674 

0.524 

0.388 

0.280 

0.200 

0.142 

Tg (IN s) 

0.187 

0.187 

0.240 

0.347 

0.560 

0.960 

0.813 

500 MHz 

AT (IN K) 

0.418 

0.301 

0.234 

0.173 

0.125 

1 0.0896 

^0.0637 

Tg (IN s) 

0.187 

0.213 

0.267 

0.347 

0.560 

0.987 

1.82 

2 GHz 

AT (IN K) 

0.209 

0.151 

0.117 

0.0867 

0.0626 

0.0448 

0.0318 

Tg (IN s) 

0.313 

0.240 

0.267 

0.373 

0.587 

0.986 

1.80 


In order to compare the sensitivities listed in table 4-2 to the 
corresponding values for the analog, system, the measurement time of the 
digital process must be defined. The step response of the control loop 
is determined from equation ( 4 - 7 ) to be 

y(nxg) = C(Tb - Ta)(1 - (1^.23) 

where is the loop sample rate (x^ = ^ . The average based on 

N samples in the post-loop processor is 
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k-N+l 


yj^(kTo) 


= c(tb - t^) 1 - i y 


-50(k-n)x, 


where yjj(kTQ) is the output of the post-4.oop processor at kr^ 


(Tq is the sample time based on the post-loop sample rate). The 
error of the step response is determined from equation (4-2U) to 


^y^( hx^) - c 


(Tb - Ta) -50(k-n)x^ 


(l^-25) 


If the settling time is again defined as the time when the error due to 
the finite settling time is less than 10 percent of the fluctuation 
error, then 


Ay (k X ) = -^ = (T - T ) i - 50 (ks-n)xo 


(i+-26) 


where h^x^ equals the settling or measurement time for the yjj post- 
loop output. The worst conditions for the settling time occur when 
T^ = 0.0 K and when the first sample is taken just as the input step 
goes high. Under these conditions, the settling time is 


_ 1 ^ - -1 0 NAT(N)/l - e^^^° 

s^ ) ^s o ’5o 3080 ( 50 XqN 

\1 0 


( 4 - 27 ) 


where AT(N) is the sensitivity of the y.. post-loop output and x (N) 


is the settling tipie of the output, 'Jhe settling times for each of 
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the post-loop processor outputs have been calculated using equa- 
tion {k-2'j) and are listed with the appropriate sensitivities in 
table k~2, A graph of required settling time versus sensitivity for 
the digital processing radiometer is presented in figure 4-2. 

Of interest also is the level of the Dicke ripple. While the odd 
harmonics of the Dicke switch frequency are eliminated by the accumu- 
lator as described in Chapter III, even harmonics will be generated due 
to the noise injection pulses. Consider the block diagram of the radiom- 
eter in figure 2-5, and the corresponding error voltage waveform at the 
output of the correlator Illustrated in figure 4-3. pue to the 
linearity of the Fourier operator, the even harmonics of this asym- 
metric waveform can be considered to originate from a pulse train at 
the Dicke frequency, with duty cycle where d = ^ is defined as 
the ratio of the time of the injected noise to one-half the Dicke 
period (see figure 4-3). The magnitudes Vjj of the sine-cosine Fourier 
expansion of this pulse train can be expressed as 


% ~ ^ ( f) I 



(4-28) 


The accTirnulator algorithm eliminates the odd harmonics, the notch filter 
eliminates the 600 Hz component, and the loop response attenuates the 
higher order even harmonics. Thus, the most significant component of 
the Dicke ripple is the fourth harmonic 


T 

^4 0 [sin (2TTd)J 


(4-29) 
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The desired signal, i.e., the do output level of the accumulator, is 
given hy 



(it-30) 


The ripple factor p expresses the output peak ripple as a percentage 
of the dc output. 


P =-x 

o 


p = 


Tj 

Kj^KgloCjajij,) I sin (2Trd) x loo^ 

BA 


(4-31) 


Substituting the appropriate value for the loop gain equation (4-31) 
becomes 


P 


9.78 X 10-4 _ 

2Trbo8 - T^) 



(4-32) 


This ripple factor versus input temperature is plotted in figure 4-4. 
The Dioke ripple level within the digital radiometep is far below the 
ripple in the analog radipmeter. In addition to the reduction in the 
level of the Dicke ripple, the digital signal processor is also more 
efficient than the analog unit in that it can form a measurement of the 
same sensitivity as the analog radiometer in less time. In the follow- 
ing chapter a more detailed comparison between the analog and digital 
signal processors will be made. 






B. Experimental 

A prototype version of the digital signal processor has been imple- 
mented on an Intel 8086 microprocessor evaluation system. The digital 
processor is shown in the user breadboard area of the evaluation system 
illustrated in figure 4-5. Only a small fraction of the support chips 
and memory available in the evaluation syston were actually used by the 
digital processor. However, the monitor program and thorough documenta- 
tion available for the evaluation system made it ideal for the first 
realization attempt. The prototype digital processor was then inter- 
connected with the existing radiometer so that it replaced the existing 
analog estimation and averaging sections (see figure 4-6). The digital 
processor was tested for proper operation as a part of the radiometer, 

The output of the video amplifier for the radiometer with the analog 
processor is shown in figure 4-T(a). The Dicke switching process can 
easily be seen. The noi$e injection pulses can also be seen during the 
antenna portion of the Dicke switch qycle. The output of the video 
amplifier with the digital processor inserted is shown in figure 4-T(b). 
The noise injection pulse can clearly be seen in the first part of the 
antenna portion of the Dicke switch cycle. The step response of the 
analog radiometer [lO] is shown in figure 4-8(a). Note the Dicke ripple 
seen in the response. The step response of the digital processing radiom- 
eter is shown in figure 4-8(b). From figure 4-8 it can be seen that the 
Dicke ripple has been reduced significantly. 
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(b) 

Figure i+-7 . Video amplifier output of analog and 
digital processing radiometers. 
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The radiometer has been installed on an aircraft, and data were 
recorded in flight with both analog and digital signal processors. 
Results of this test flight are very encouraging. Analyses performed 
on these test data thus far indicate that the digital processing 
radiometer is operating as predicted. 
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CHAPTER V 

CONCLUSIONS ^ 

The p-ur pose of this research was to develop an improved power estima- 
tion strategy for a Dicke switching noise injection radiometer. It was 
shown that the Dicke ripple inherent in the Dicke switching process 
within the analog radiometer could be eliminated using sample data 
techniques. Elimination of the ripple allowed the use of a faster 
responding control loop. The resulting decreased settling time leads to 
an increased number of uncorrelated samples at the loop output, thus 
allowing most of the averaging to be performed in a post-loop processor. 

The flexibility and efficiency of this post-loop averaging result in an 
improved estimation strategy. A comparison of the radiometer using the 
proposed processing technique and the original analog radiometer is pre- 
sented below in support of the improved performance. 

A. Dicke Ripple 

An expression for the Dicke ripple level as a percentage of the 
desired dc signal was determined for both the digital and analog pro- 
cessors in previous chapters. It was determined that this ratio was a 
function of antenna temperatiire. A graph of the ratio of the ripple 
level to dc level (or ripple factor) versus antenna temperature is shown ^ 

in figure 5-1. In this graph, the ripple factors for the digital and 

« 

analog processes are shown on the same scale. The reduction in the 
ripple level due to the elimination of the odd harmonics of the Dicke 
frequency by the accumulator can clearly be seen. 




Figure 5-1. Comparison of ripple factor 
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This reduction in Dicke ripple is significant since it permits the 
control loop to have reduced settling time. This fast settling time is 
available since the loop is no longer required to attenuate the Dicke 
ripple. 

B. Sensitivity 

As explained in Chapter II longer Integration times produce a more 
sensitive measurement of a larger area (less spatial resolution). It i^ 
then desired to perform the measurement to the required sensitivity in 
the shortest possible time. Thus, in order to compare the sensitivity 
of the analog and digital radiometers, the time required to form the 
measurement must be considered. The Y 2 post-loop output, or the 
average of two loop output samples, has approximately the same theo- 
retical sensitivity as the original analog radiometer. The times 
required to form the measurement of these sensitivities for the avail- 
able input bandvidths are shown for the analog and the J 2 post-loop 
output of the digital processing radiometer in figure 5-2. From fig- 
ure 5*“2, it can be clearly seen that the digital processing radiometer 
is superior to the analog unit in that it can provide a measurement with 
the same sensitivity as the analog processor in less time. This results 
in better spatial resolution. 

C. Spatial Resolution 

For a specified spatial resolution the total error in the estimate 
of the brightness temperature is made up of fluctuation error and 
^^antenna smear” error. Longer integration times will reduce the 
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Figure 5-2 • Comparison of req^uired measurement time 




fluctuation error AT. However, longer integration times will degrade 
the spatial resolution, due to the measurement platform movement, thus 
increasing the "antenna smear" error given in equation ( 2 - 15 ). Examina- 
tion of equation (2-15) reveals that the "antenna smear" error is a 
function of the brightness temperature. Referring to equation (2-15) 
it can be seen that if the radiometer is moving over a scene in which 
the brightness temperature changes slowly, it is possible to increase 
the integration time with little increase in the "antenna smear" error. 
However, as the radiometer approaches a point where the brightness 
temperature changes rapidly (such as an ice/water boundary), the inte- 
gration time must be reduced to limit the "antenna smear" error. 

An advantage of the digital processing radiometer is that the 
measurement is made to several sensitivities simultaneously, and these 
estimates are available simultaneously at the various post-loop outputs. 
Thus, an operator can determine which post-loop output will provide the 
best estimate of the brightness temperature for a particular type of 
scene and specified spatial resolution. More specifically, for a 
desired spatial resolution, there is an optimum integration time which 
will minimize the total adverse effects due to fluctuation errors and 
the "antenna smear" error. This decision must be made by the operator 
with the present system. However, the antenna footprint size and air- 
craft speed may be such that it could be possible to approximate the 
infinite sum of equation (2-15) by the first few terms of the series. 

This would result in an approximation of the "antenna smear" error for 
a particular footprint size, aircraft speed, and scene. Then, using the 
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radiometer's estimate of the scene, a continual approximation of the 
"antenna smear" error could be formed. If this is possible, the micro- 
processor program! could determine the optimum integration time for a 
particular required spatial resolution. Thus, the radiometer could con- 
tinually adjust the integration time as the scene varies by selecting 
the appropriate post-loop output such that the output is the best esti- 
mate of the scene for the required spatial resolution. This varying 
integration time scheme is presented here to demonstrate the desirability 
of the added flexibility obtained using the digital processor. Although 
no further investigation was attempted at this time, it is suggested 
as a topic for future research. 

D. Summary 

The digital signal processing radiometer presented here is superior 
to the original analog system in several respects. The Dicke ripple due 
to the modulation and demodulation processes within the radiometer is 
substantially reduced by the digital signal processor, thus enabling the 
use of the post-loop processing technique. The time required to form a 
measurement to a specified sensitivity is shorter for the digital pro- 
cessing radiometer. This reduction in required measurement time is 
actually due to the more optimum loop design and post-loop averaging, 
and is not specifically a result of the digital techniques used. 

However, the subsequent reduction in measurement time is significant 
since it results in a more sensitive radiometer. The added flexibility 
obtained by use of the digital techniques may also lead to improved 



radiometer performance. Thus, the objective of this investigation, i.e 
the development of a power estimation strategy demonstrating perforraanc 
superior to the original design has been accomplished. 
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APPENDIX A 


OVERVIEW OF THE DICKE RADIOMETER 


A radiometer is a precise instriunent which measures electromagnetic 
energy radiated from a surface. Any object will emit a particular 
broadband spectrijm of electromagnetic energy depending on its tempera- 
ture and physical characteristics such as emissivity. Max Planck 
established that the emission spectrum of a hypothetical blackbody is 
completely governed by its absolute physical temperature [j]. The power 
spectrum of a blackbody with an absolute temperature T can be 
expressed by 


P 


' [exp(hf/kT) - ^ Af dJ2 


(A-1) 


where 

P = power radiated into solid angle dJ2 in bandwidth Af, W 
h = Planck's constant = 6.63 x 10~^^ J*s 
f = frequency, Hz 
c = velocity of light, m/s 
k = Boltzmann's constant = 1.38 x lo”^^ J/K 
T = absolute physical temperature, K 
In the microwave frequency region, for most practical situations, the 
term hf/kT is sufficiently small that the exponential function can be 
approximated by the first term of its Taylor expansion; i.e., 
exp(hf/kT) - 1 + (hf/kT). The power radiated by the blackbody can 


then be expressed as 



87 


P = Af 6.Q 


(A-2) 


where X = c/f. 

An ideal "lossless" linearly polarized antenna with receiving 
cross section A(ft) viewing a blackbody receives a power given by 


P =1 

R 2 


i 


2kT Af 


A(fl) dfl 


(A-3) 


where the factor of one-half is due to the fact that the antenna receives 
energy in only one plane of polarization [ 7 ]. 

It can be shown [2] that 


2 

r A{Q) df2 = f G(fl) dJ2 

Jq 


4Tr 


= 


ikJl) 


U-h) 


where G(fil) is the absolute antenna gain function. Combining equa- 
tions (A-3) and (A-4) results in 


Pj^ = kT Af 


Rearranging equation (A-5), 


(A-5) 


T = 


R 


k Af 


(A- 6 ) 


The quantity T is called the antenna temperature, which will subse- 
quently be denoted as T^. The value T^ corresponds to the temperature 
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that a hlackbody must have at the antenna terminals to produce the 
received power The rationale for the common practice of referring 

to a noise power distribution as an equivalent noise temperature is 
illustrated by equation (A-6). 

In practice, the object of interest is not a blackbody. Thus, 
the received power is not solely a function of the physical temperature, 
but is also a function of the emissivity, reflected radiation, and 
atmospheric path. This investigation only considered the problem of 
accurately measuring the antenna output power or antenna temperature 
T^. The problem of determining properties of the object from the 
measurement of T^^ was not considered. 

The output signal from the antenna is actually low level gaussian 
noise with zero mean. The radiometer must estimate the variance or 
noise power of this signal. Measuring the variance of this signal pre- 
sents several interesting problems. One problem is that the actual 
level of the received signal is extremely low and is often smaller than 
the receiver noise. The antenna temperatures of interest range from 

100 to 270 K, corresponding to a received power density of I.38 x lO"^^ 
-21 

to 3.73 ^ 10 W/Hz. In comparison, the noise generated within the 

radiometer referred to the input is approximately 627 K, corresponding 

—21 

to a receiver noise power density of 8.65 ^ 10 W/Hz. Thus, the 
desired signal appears to be obscured by the higher level receiver 
noise. Also, due to the extremely low level of these signals, sub- 
stantial amplification will be required before any attempt can be made 
to measure the level of the received power. The amplifiers used to 
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provide this gain must be very stable since changes in the gain will 
cause changes in the measured signal power. Thus, in order to measure 
the power of the received signal, the receiver noise and the gain of 
the radiometer must be precisely known and extremely stable. Any 
uncertainty in the receiver noise power or the gain will result in an 
erroneous measurement of the received signal power. 

The Dicke switching noise injection radiometer shown in figure A-1 
is one system used to measure the power of the received signal while 
reducing errors due to the previously mentioned problems. A brief 
introduction to the basic operation of this radiometer will now be 
presented. Referring to figure A— 1, the input to the radio frecjuency 

(RF) amplifier is switched between a known noise temperature T and 

« B 

the antenna output plus feedback injected noise T. + T_. This switched 

A F 

signal passes through the RF amplifier and input filter where it is 
corrupted by the receiver noise T^. Thus, the output of the filter 
switches between two zero mean gaussian distributed noise signals, one 
with variance proportional to T^ + Tj^, and the other with variance 
proportional to T. + T„ + T„. 

A i* R 

The output of the filter is applied to the square-lav detector. 

If the input to the square-lav detector is a random process x(t), then 
the output is given by 

y(t) = x^(t) 


(A-7) 













The mean of y(t) is given by 


My(t) = E{x^(t)} 

yy(t) = [E{x(t)}]^ + a^(t) (A-8) 

u (t) = y^(t) + a^(t) 

^ X A. 

Since ~ equation (A-8) can “be reduced to 

yy(t) = 0^(t) . (A-9) 

Thus, if the input to the square-lav detector has zero mean, the mean 
value of the output of the square-lav detector equals the variance or 
pover of the input. The input to the square-lav detector is svitched 
hetveen tvo zero mean random variables, one vith variance proportional 
to Tg + and the other vith variance proportional to T^ + T^ + T^. 
Thus, from equation (A-P), the output of the square-lav detector 
svitches hetveen tvo random variables, one vith mean proportional to 
Tg + T^, and the other vith mean proportional to T^ + T^ + T^^. The 
output of the square-lav detector is svitched in synchronization vith 
the input svitch such that the signal is alternately connected directly 
to the integrator and then inverted before being connected to the inte- 
grator. When both svitches are up, the integrator ^^sees” a random 
variable vith mean proportional to T^ + T^j, + T^^. When both svitches 
are dovn, the integrator "sees” a random variable vith mean proportional 
to ""(Tg estimation and averaging section then estimates the 



mean (i.e., dc value) of this random variable. The mean of the signal 
qver a full switching cycle is given by 


92 


We = C[(T^ . Tp + T^) - (Tg . T^)] 


We = 0(T^ + Tj, - Tj) 


(A-10) 


where C is a system constant. It can be seen from equation (A-IO) 
that the receiver noise does not contribute to the net mean value of the 
input to the estimation and averaging section. Thus, the Dicke switch- 
ing process effectively eliminates the need to know the value of T^^, 
as well as the uncertainty in the measurement due to fluctuations 
of Tj,. 


The output of the estimation and averaging section is then used to 
control the injected noise. Since the noise tanperatures of the refer- 
ence source Tg and injected noise T^, are known, an estimate of T^ 
can be made from the result of equation (A-IO). This estimate is used 

to adjust the injected noise power such that T„ = T_ - T. . The mean 

r n A 

value of the input to the estimation and averaging section (or loop 
error signal) as determined from equation (A-IO) then becomes 


Pe = C[(T^ + Tp . Tp) - (Tj + T^)] 

Pe = c[(T^ . (Tj - T^) . Tp) - (Tp . T^)] 

Pe = c[(Tb + Tp) - (Tp . Tp)] 


(A-11) 


0. 
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Thus, the mean vaJ-ue of the random variable input to the estimation and 
averaging section is zero, and the control loop is locked. This feed- 
back concept will automatically compensate for fluctuations in the 
forward loop gain since such changes will not affect the final value of 
the loop. Thus, the Dicke noise injection radiometer has eliminated the 
effect of the receiver noise power and has reduced errors caused by 
fluctuations in the forward loop receiver gain. 





APPENDIX B 


IMPLICATIONS OF TIME AVERAGING A NONERGODIC PROCESS 
In order to determine what the time average of the output of the 
square-law detector implies about the ensemble average, consider the 
random process y(t) at the output of the square-law detector. Let 
Viy(t) represent the time average of y(t), and let 

y(t) = 1iy(t) + z(t) (B-lj 


where z(t) is a random process with zero mean, and y^(t) is the mean 
value of y(t) corresponding to the brightness temperature. The time 
average of the output of the square- law detector is then determined to 
be 


A 1 

r t+Tl2 


yy(t) = - 

1 

y(C) dc 

“ L 

4 

■'t-T/2 



r t+x/2 

[yy(c) + z(c)] d? 

P/t) 

't-x/2 


(B-2) 


If this time average is an unbiased estimate of y (t), i.e., the scene 
brightness temperature, then E|yy(t)| should equal the brightness 
temperature B(t). Taking the ensemble average of both sides of equa- 
tion (B-2) results in 
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* 


A 


r, f t+x/2 

[Uy( 0 ..( 0 ] 



®| Sy (*)| = 7 / 


t+T/2 


E)iiy(t)| 

K|V*)1 


= i f** 

^ •'t-T 

--T 

^ h-T 


t-x/2 
t+x/2 
/2 
t+x/2 
/2 


E)yy(c) I + e|z(c) { dC 

Ejyy(?)( dC 

B(c) dC # B(t) 


(B-3) 


From equation (B-3) it is clear that the time average is not an unbiased 
estimator of the brightness temperature. However, a closer look reveals 
that the time average is an unbiased estimate of the average brightness 
temperature of the scene covered by the antenna footprint in T seconds. 

If the time average is used as an estimator of the actual bright- 
ness temperature within the antenna footprint, a bias error will result. 
This bias error can be obtained by expanding the right side of equa- 
tion (B-3) in a Taylor series. Performing this expansion about ^ = t 
results in 


i 





A 1 r / \ 2 

E|yy(t)f = - I B(t) + (? - t)B'(t) + B"(t) 

•'t-r/2 

+ B'"(t) + . . . dc 


.t+T/2 


Ejy (t)[ = B(t) + B'(t) if* (? - t) d? 


/\ 2 ^ 

E|y (t) [ = B(t) + B"(t) |r- + B""(t) 


(C - t)^ 

't-x/2 

t£ 

24 


dC + . . 


1920 


(2n) 2n 


E|y (t) [ = B(t) + y B(i)_L 

j^Ti (2n)!(2n + 


1 )( 2 ^'^) 


(B-4) 


where B'(t) and represent the first and 2nth derivatives 

of B(t) with respect to t respectively. The infinite sum of 
equation (B-4) is the resulting bias error if the time average is used 
as an estimate of B(t). This bias error is due to the fact that the 
process is nonergodic. 
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